Signature of stripe pinning in optical conductivity 
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The response of charge stripes to an external electric field applied perpendicular to the stripe 
direction is studied within a diagrammatic approach for both weak and strong pinning by random 
impurities. The sound-like mode of the stripes described as elastic strings moves to finite frequency 
due to impurity pinning. By calculating the optical conductivity we determine this characteristic 
energy scale for both a single stripe and an array of interacting stripes. The results explain the 
anomalous far-infrared peak observed recently in optical-conductivity measurements on cuprates. 
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Recent optical conductivity measurements in high-X" c 
cuprates have revealed the existence of a strong peak 
in the far-infrared regime, which has been attributed to 
the presence of charge stripes Q, 0] . The typical energy 
of this peak is much lower than the mid-infrared band 
detected in previous experiments Q and reproduced by 
numerical simulations of the Hubbard model 4]. A sim- 
ilar feature reported for slightly overdoped B^S^CuOe 
[f| was explained in terms of a charge-ordering instability 
of the Fermi liquid @]. However, these arguments cannot 
be extended to the underdoped regime, where the Fermi- 
liquid description breaks down and the charge degrees of 
freedom are known to order [3. 

In this paper we propose that this anomalous peak may 
originate from the pinning of the zero-energy phonon-like 
mode associated with transverse fluctuations of stripes, 
which shifts to finite energies in the presence of impuri- 
ties. We calculate the optical conductivity of striped sys- 
tems using a diagrammatic approach, which allows us to 
determine the peak frequency for both weak and strong 
impurity pinning. In the single-stripe case, we find a peak 
at a frequency v given by the ratio between the stripe 
velocity v and a characteristic impurity length scale A, 
which is determined by the strength of the pinning. For 
a stripe array, a second energy scale lujj, proportional to 
the strength U of the stripe-stripe interaction, becomes 
relevant. In the weak-pinning regime, which is appropri- 
ate for describing cuprates, a second peak appears around 
u>u- When u>u ~ v, a resonance arises and the main 
peak splits. Our calculations yield a good estimate for 
the peak frequency measured in La cuprates 0, 0, @ and 
clarify some peculiar features observed in experiments, 
such as the peak-sharpening and splitting which 
occur when stripe-stripe interactions are important. 

We begin by considering a single line of holes (stripe) 
embedded in a two-dimensional i x L antiferromagnetic 
background with lattice constant a. The transversal di- 
mension I correspond to the interstripe distance when 
the stripe array is explicitly considered (see below). The 
stripe is oriented along the y-axis and acts as a domain- 
wall for the underlying magnetization. The holes are as- 
sumed to move in the transverse x direction and their 



dynamics is governed by the t — J model. In the long- 
wavelength limit a wave-like action describes the trans- 
verse displacement u(y) of the stripe with respect to the 
equilibrium position 9]. The stripe velocity v = ^/jia 
is related to the kinetic energy of the holes (t) and the 
staggered magnetization (J) of the surrounding regions. 
In terms of the dimensionless displacement field ip = u/a, 
the action reads 



S =2^/ drdy [C Q ty] 



(i) 



where Co[ip] = (dt/j/dr) 2 + v 2 (dip/dy) 2 and V[ip] is the 
pinning potential. In analogy with the standard prob- 
lem of pinning of elastic manifolds 0, 0, 0] , the inter- 
action of the (striped) density of carriers p(r) with the 
short-range disorder potential U(r) = Uq J2i <5( r ~ r i) can 
be written as / drU(r)p(r) ~ poUo J2i cos[27r(u — Xi)/£], 
where r; = (x,,^) denotes the position of the defects, 
and po is the average stripe density 0]. We consider 
that only defects localized near the stripe effectively act 
to pin it, so that the cosine term can be expanded around 
the minimum to yeld, within the quadratic approxima- 
tion, V[ijj] = V atJ2 i ('4 ; (y)-Pt) 2 ^(y-yi)- The impurities 
pin the local displacement ip(yi) at a value Pi, with pos- 
itive strength Vb, and we assume % = ks = 1- 

In the absence of dynamical fluctuations, the stripe 
accommodates the impurity potential by adopting a con- 
figuration ipo(y), which is a solution of the static equation 
of motion. The full, time-dependent solution ip(y,r) is 
described as a fluctuation around the equilibrium config- 
uration -0o (y) , i.e., ip(y,r) = ip (y)+<P(y,T~). Introducing 
the Fourier transform <j)(q,u> m ) with respect to the mo- 
mentum q and the Matsubara frequency u> m = 2-KmT ', 
the reciprocal-space action for the cf> field becomes 



S = - V o 1 ( ( li ul m )\(j)(q,Lj m )\ 2 



To 



(2) 



9,9' 1^7. 



where S(q) = (1/L) £i 



is the form factor, and 



V (q,uj m ) = at(u)m + v q ) is the bare Green function. 
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An electric field applied in the x-direction, perpendic- 
ular to the stripe, generates a current J = J dye(d4>/dt). 
The linear charge density of the stripe is denoted by el a . 
From the Kubo formula for the optical conductivity |13| , 



a(uj) = e iujV(q = 0, iui r , 



id) and by analytic 



continuation in the lower half-plane of the dressed Green 
function T>(q,Lu m ) one finds 



%ta(uj) = —e 2 ui^sT>(q — 0, iuj m — > co — id). 



(3) 



In the absence of impurities, T> coincides with the bare 
Green function Vq. Its replacement in Eq. yields 
the expected response of a massless, sound-like mode, 
3?cr(w) = e 2 atd(uj). The dressed Green function T> may 
be evaluated from the Dyson equation, T>^ 1 (q,uj m ) — 
Vq (q,uj m ) — £ = (uj^ + v 2 q 2 — Ft) /(at), where we have 
performed a rescaling S — Y/a such that Y has dimen- 
sions of energy. The self-energy Y is evaluated by av- 
eraging the product of the several S(q) factors, which 
appear in each term of the perturbative series for T, over 
the random impurity positions |l3j . Crossing diagrams, 
arising from scattering processes by different impurities, 
correspond to increasing powers of the impurity density 
n$. If r is evaluated up to order Vq 2 , crossing terms do 
not appear (Born approximation) and one obtains 



Y(uj m ) = -~riia + ( -y J n$ay ^D (g,q; m ) 



(4) 



where Vq was used for calculating the second-order cor- 
rection. On defining the frequency ujq — mv = maVJt, 
the dimensionless quantities uj = uj/luq, G — Yt/ui 2 ,, and 
the parameter a = Vo/(2maJ), Eq. J3J takes the form 



G = 



2Cj 



(5) 



According to Eq. JjJJ the optical conductivity is given by 

CjG" 



(u 2 + G') 2 + G" 2 ' 



(6) 



where <tq = e 2 at/u>o and G = G'+iG". Note that as soon 
as a real part is generated in G, the delta function in the 
optical conductivity moves to a finite frequency and its 
amplitude is determined by G". Using (JjJJ, we find 



5Jct(w) = (To 



2(ujvf 



A(lo/v) 2 [(uj/v) 2 - l} 2 + a : 



(7) 



with v — y/auiQ. We note that 3ta(ui) — > both as u> — > 
and as lu — > oo, and has a peak at lu ~ ^, as shown by the 
solid line in Fig. la. The parameter a, which controls the 
perturbative expansion for G, depends both on the ratio 
Vq/ J, i.e. the competition between impurity potential 
and elastic energy, and on the number of impurities per 
site, riia. Thus the weak-pinning limit, where only first 
powers of a need be considered (as in Eq. ©), is realized 



for a large number (ma ~ 1) of weak (Vb/J <C 1) impu- 
rity centers. In this case, Eq. Q shows that the peak 
frequency scales according to v ~ ^/Votma, moving to- 
ward zero energy as the strength of the pinning potential 
decreases. In the evaluation of Eq. (0J, we used the bare 
Green function T> . The divergence of the self-energy (JSJ 
at small lo indicates the break down of this crude approx- 
imation at low frequencies. Inclusion of the next-order 
term in the perturbative expansion shows that higher- 
order diagrams become important for to < i.e. 
at frequencies much smaller than v. It has been argued 
[ll| that a better result for G may be achieved by using 
the dressed Green function in Eq. JSJ) , which leads to the 
self-consistent equation G = —a+ (a 2 /2)(—uj 2 — G)~ 1//2 . 
However, the resulting expression for a (to) shows a spu- 
rious gap below to w 0.7i>, which is thought to originate 
from the fact that the self-consistent Born approximation 
selects arbitrarily a subset of diagrams of higher order in 
a without including crossing terms, whose contribution 
cannot be excluded if ma is large. Because we are inter- 
ested primarily in the evaluation of the peak frequency z/, 
and not on details of the low-frequency behavior of o(us) 
[T^ |. we restrict our analysis to the non-self-consistent 
Born approximation, which provides a reliable estimate 
both of the maximum of <j(u>) and of its scaling with a. 

A different approach must be considered in the strong- 
coupling limit a>l, which corresponds to a small num- 
ber (ma, <C 1) of strong (Vq/J ^> 1) pinning centers. In 
this case one might expect that the entire perturbative 
series in a should be resummed, but there exists a cri- 
terion for selecting the relevant higher-order diagrams. 
For ma <§; 1, it is possible to select the subset of dia- 
grams of first order in m after averaging over the impurity 
distribution (T- matrix approximation) [l3| . By explicit 
summation of the resulting series, using the full Green 
function V while evaluating each diagram, one obtains a 
self-consistent equation for G. In the limit a — > oo, one 
finds the analytical solution 5Rcr(w) = 2an Vuj 2 — 1/oj 3 for 
u> > 1 and 3?cr(w) = for w < 1 12]. Thus, Sa(w) 
displays a gap at ujq and a maximum at a frequency 10 
slightly above this value. When the solution is evaluated 
at arbitrary a 1, one finds essentially the same fea- 
tures. In this case the requirement of full self-consistency 
is crucial, and the gap feature arises from the strong 
pinning of the sound mode. The difference between the 
weak- and strong-pinning results can be understood by 
considering that in both cases the peak frequency cor- 
responds to the sound mode of a free string, which is 
now trapped on a characteristic length scale A. When 
few strong impurities are present, A is given by the av- 
erage distance between impurities 1/m and the peak ap- 
pears at the frequency luq. However, when the number 
of impurities increases and Vq decreases, the elastic en- 
ergy cost of accommodating the string on a length scale 
l/n$ is too high, and ip is pinned on the Larking length 
A;„c 0, E|- Dimensional estimates based on Eq. (JTJ 
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yield A; oc ~ v/y/Voriiat, which indeed corresponds to the 
characteristic frequency v in the weak-pinning regime. 
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FIG. 1: Left: real part of the optical conductivity, in units of 
(To, as a function of u/v evaluated in the Born limit (a — 0.1). 
When interaction between stripes is not relevant (wp < v, 
Fig. la) one finds a single peak (the second one for the stripe 
array is almost suppressed), as observed in LSCO at x = 0.05, 
Fig. lb Q. When the resonance condition luu = v is satisfied, 
the peak at u — v splits into two equally-weighted peaks (Fig. 
lc), as observed experimentally at x = 0.07, Fig. Id |3- By 
increasing the inter-stripe interaction (uu = 1.5f, Fig. le) 
the peak sharpens and moves to slightly lower frequencies, 
as observed experimentally by comparing Nd-free with Nd- 
doped LSCO compounds, Fig. If For comparison with our 
zero-temperature calculations, we display the experimental 
data at the lowest measured temperatures above T c . Notice 
the logarithmic scale on the right panels. 

Until now we have concentrated on the transverse pin- 
ning of a single stripe. However, a more realistic de- 
scription of cuprates requires the study of an array of 
stripes. By considering a harmonic coupling strength U 
between neighboring stripes, which are labeled by an in- 
dex n — 1, . . . N, we obtain the action |l6( 



S = ^ £ / dTdy + Ut ^ n ~ ^" 



Here tp„ — i/) n (y,T) and the fluctuating held <p„{y, r) = 
(f>{x — n£,y,r) is a function of the discrete x values. 
By repeating the same steps as before, we observe that 
Eq. J2J is still valid, on taking q — > q = (q x ,q y ), 
S{q) = (LN)- 1 Y Jm e tq:cnl+lqvy * ", where yf is the impu- 
rity coordinate at the n-th stripe, and T>q becomes 

at 



£>o(q, u m ) 



v 2 ql 



AUt sin 2 (qJ/2) 



(8) 



The system of interacting stripes displays an anisotropic 
sound-like mode at long wavelengths, with an effective 
elastic coefficient U/a in the x direction, and an optical 
mode at q x £ = n with energy u>u — 2yUt. 

The general procedure used above to address the effect 
of disorder is extended readily to the case of interacting 
stripes. We assume that the average density of impurities 
does not depend on the stripe index n, i.e. n™ = n^. In 
the Born limit, the main difference with respect to the 
single-stripe case arises from the additional integration in 
q x of the Green function T>q in Eq. J3J. By introducing 
the rescaled variables and Q = lu/lou, one obtains 



7T UU Jo \/ 



1 



sm 2 x 



(9) 



-i) 2 



When u> > 1 the second term in Eq. @ contributes only 
to G" . In the limit of vanishing stripe-stripe interaction 
U — * 0, or u> — * oo, it reproduces the single-stripe re- 
sult Eq. |JSJ. For finite U, a second energy scale uou is 
relevant for determining the optical conductivity, which 
is defined by Eq. ((3J) with q replaced by q. The new 
feature is the divergence of G" at uju, which leads to 
a second peak around this frequency, where $ta(u)) van- 
ishes. The vanishing of the conductivity at lo\j in the 
weak-pinning regime should indeed be expected: at this 
energy, neighboring stripes fluctuate in antiphase, and 
because the total current arises from the sum of the cur- 
rents of each stripe, the net current in such a case is zero. 
At uju <C v or luu v the peak at u)jj is almost com- 
pletely suppressed by the overall decrease of dt<r(uj) at 
small and high frequencies, see dashed lines in Figs, la 
and le. However, when ujjj ~ v the peak at lu ~ v splits 
(Fig. lc), signaling a resonance between the characteris- 
tic mode of the stripe array at the zone boundary and the 
characteristic frequency of sound-like oscillation of each 
stripe at the Larkin length. Dimensional estimates indi- 
cate that this resonance arises when the average stripe 
separation coincides with the Larkin length. Moreover, 
for the stripe array we observe a sharpening of the peak 
for lo\j > v due to the flattening of G" around v (dashed 
line in Fig. le). The softening of the peak frequency for 
uju > v scales according to oPv/uu, so it never moves 
far away from v. 

In the strong-pinning regime the qualitative behavior 
of 3ta(uj) for the stripe array is similar to that of the 
single stripe, and no resonance is observed. Because the 
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interaction with the impurities is the predominant effect, 
the behavior of G is almost entirely dominated by the 
pinning of the sound mode at ojq. More details about the 
strong-pinning regime will be presented elsewhere [l7| . 

We next compare our results with the available exper- 
imental data. The velocity of the unpinned sound mode 
may be estimated by using standard values t ~ J ~ 0.1 
eV for cuprates. The other relevant quantity is n,. In 
these materials, the chemical-substituted dopants, whose 
two-dimensional density is x/a 2 (x being the doping), si- 
multaneously provide charge carriers in the plane and act 
as pinning centers. As we explained at the beginning, we 
consider that on average the impurities within a range 
£ contribute to pin each stripe, i.e. n,j = x£/a 2 . Neu- 
tron diffraction measurements have shown that in the 
underdoped regime (0.05 < x < 0.12) of La2- x Sr x Cu04 
(LSCO) t = a/2x 7], which yields n ia = 1/2 and 
uio ~ 0.05 eV. Because n^a ~ 1 corresponds to the weak- 
pinning regime, we expect that the peak frequency is 
given by v — \/au>o and thus is reduced with respect to 
the estimate of loq, in agreement with experimental data, 
see Fig. lb, Id, If As shown in Fig. lb-lf, the 

peak frequency softens as x increases. Even though we 
also expect a general decreasing of v due to a screening 
effect on Vo, a more quantitative comparison with experi- 
mental data would require the inclusion of lattice effects, 
which is beyond the scope of this paper. 

Another interesting feature can be explained by 
our calculations: with increasing doping concentration, 
stripe-stripe interactions become more important and 
one should consider the stripe array. Our results predict 
that the peak splits when the average stripe separation 
£ is of the order of the Larkin length, see Fig. lc. This 
splitting has been observed: for x = 0.05 doped LSCO, 
a single peak appears (Fig. lb) whereas for x — 0.07 
the peak splits (Fig. Id) [2|. Although the value of a 
as a function of doping x is in general not known, we 
may nevertheless estimate the peak frequency where the 
splitting arises as v — vjl — 0.014 eV = 113 cm -1 , in 
excellent agreement with the measured value (Fig. Id) 
■ Our calculations also allow us to understand the re- 
sults obtained by Dumm et al. 1] in Nd-doped LSCO 
(Fig. If): the broad peak observed in LSCO moves 
to slightly lower frequency and becomes sharper in the 
Nd-doped compound, where the stripe array is ordered 
0- We have shown that when inter-stripe interactions 
increase, the peak becomes slightly softer, narrower and 
more intense (Fig. le), as observed experimentally. Note 
that in Ref. |lj the feature seen in the Nd-free sample 
in Fig. If was interpreted as an anomalous large Drudc 
peak. However, we claim that a stripe signature is still 
present in these data, even though a partial overlap with 
the single-particle response is possible. Indeed, resistiv- 
ity data from Ref. give tr(0) ~ 8- 10 3 fi -1 cm -1 , which 
is much larger than the value obtained by extrapolating 
the supposed Drude peak. Moreover, a peak at nearly 



the same frequency is also seen at x = 0.12 in Ref. |2|. 

In addition, our results provide the correct order of 
magnitude for the peak height. The one-dimensional 
optical-conductivity unit Co = <j}f converts into a three- 
dimensional one as cr 3 ^ = a\ d jld = (2e 2 /%)(t/uo)(x/d), 
where d is the interlayer spacing. At x = 0.07, using 
d ~ 6 A, we estimate the maximum of a(u>) to be of 
order 10 3 f2 _1 cm _1 , as experimentally observed Q. 

Note that longitudinal transport, which is possible for 
half-filled stripes as in the cuprates, was neglected in our 
approach. However, the experimental observation of the 
same far-infrared peak in nickelates |l8j |. which exhibit 
filled stripes, supports our hypothesis that the anomalous 
peak in the far-infrared regime originates from transverse 
fluctuations. In fact, our results could be used also to 
describe the nickelates, but in this case the effective t and 
J parameters within a one-band model are not known. 

In conclusion, we calculated the contribution from 
stripe pinning to the optical conductivity in cuprates 
within a diagrammatic approach accounting for the pres- 
ence of randomly distributed impurities. We found a 
peak in the far-infrared range at a frequency which has 
the same order of magnitude as observed in the exper- 
iments. In addition, our results explain the splitting of 
the peak, as well as the narrowing and increased intensity, 
which arise when stripe-stripe interactions are relevant. 
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